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We present a general phenomenological model for superconducting qubits subject to noise pro- 
duced by two-state fluctuators whose couplings to the qubit are all roughly the same. In flux qubit 
experiments where the working point can be varied, it is possible to extract both the form of the 
noise spectrum and the number of fluctuators. We find that the noise has a broad spectrum con- 
sistent with 1/f noise and that the number of fluctuators with slow switching rates is surprisingly 
small: less than 100. If the fluctuators are interpreted as unpaired surface spins, then the size of 
their magnetic moments is surprisingly large. 

PACS numbers: 03.65.Yz,85.25.Cp, 85.25.Dq 



I. INTRODUCTION 

Superconducting qubits based on Josephson junctions 
are promising candidates for quantum information pro- 
cessing [l|, 0] . Integrated-circuit fabrication technologies 
provides a relatively straightforward route to scale up the 
number of qubits, and the qubit coherence times have 
been prolonged dramatically since the superconducting 
charge, phase and flux qubit designs were first developed 
over a decade ago @-|(|. However, detailed mechanism of 
decoherence due to the coupling of the Josephson device 
to external noise sources is still not fully understood [l[ . 

Recent experiments on superconducting qubits show 
that 1 /f flux noise is an important source of decoherence 
@4j| ■ Experiments over the years have agreed on certain 
universal characteristics of this noise: (1) it has weak de- 
pendence on a wide range of parameters such as SQUID 
loop geometry, inductance, material, etc.; (2) it has an 
approximately 1 /f noise power spectrum and the magni- 
tude ranges from 0.01 — lOO(fi^o) 2 /Hz at the frequency 
1Hz, where $o = h/2e is the magnetic flux quantum 

The origin of this low-frequency noise at milli-Kelvin 
temperature has been a puzzle for over 20 years and 
is still under active debate [l2T - fl4l ]. There are indica- 
tions that a high density of unpaired surface spins on 
the SQUIDs may be the physical causes of the noise 
P El EMI. These defect sites behave as two-state 
fluctuators that switches between their two states due to 
thermal activations and/or other interactions. 

In this paper, we present a phenomenological model of 
the fluctuators. The physical parameters of the model 
can be extracted from qubit measurements at different 
working points. Analysis of experiments 0, H| produces 
estimations of the effective magnetic moment and noise 
power spectrum density that are comparable to the ex- 
perimental findings. Our chief new result is that the num- 
ber of slow fluctuator is small, less than 100 and possible 
even of order 10. 

The paper is organized as follows. In Sec. UH we de- 
scribe and solve the model. This gives results for free 



induction decay (FID), energy relaxation (ER) and spin 
echo (SE) signals. In Sec. IIII[ we summarize our assump- 
tions for the flux qubit systems and demonstrate how to 
extract the physical parameters of the model from exper- 
imental data. In Sec. IIVI we discuss the results. 



II. NOISE MODEL 

The superconducting flux qubits consist of a supercon- 
ducting loop with three Josephson junctions Q . The two 
relevant states are the clockwise and counter-clockwise 
persistent current states in the loop and the loop is ef- 
fectively a quantum two-level system. The Hamiltonian 
of the superconducting flux qubit can be written as [l[ 



H qb - ~-o- z 



A 1 ( , . 



(1) 



where e and A are the energy difference and tunneling 
splitting (Josephson coupling) between the clockwise and 
counter-clockwise current states, h(t) is the flux noise 
in the environment and cr XiViZ are the Pauli matrices. 
The energy difference is proportional to the applied flux 
through the superconducting loop 



e = 2I P ($ext - *o/2) , 



(2) 



where I p is the persistent current and <f> cxt is the ex- 
ternally applied magnetic flux in the loop. When <£> e xt 
is half a flux quantum, the two current states are de- 
generate in energy. The flux noise h(t) is described by 
a time-dependent classical field. The eigenenergy of the 
qubit is thus 



B = \Te 



A 2 



(3) 



The angle 9 = tan _1 (A/e) is related to the working point 
of the device: 9 = ir/2 is the optimal point and 9 = is 
the pure dephasing point. 

The flux noise is induced by an ensemble of fluctua- 
tors, all fluctuating independently, giving rise to random 



2 



telegraph noise (RTN). Assuming a total number of K 
fluctuators, the Hamiltonian can be written in the fol- 
lowing form after a basis transformation 



1 1 K 

H = ~-B a z - 2 X! Sk ^>9k 

k=l 



(4) 



Here we redefine the z-axis to be the eigenenergy axis, gk 
is the coupling constant of the fc'th fluctuator. Note all 
fluctuators have the same 9 value since flux fluctuation is 
along the e direction. Sk (t) is the random time sequence 
due to the fc'th fluctuator and switches between the two 
values —1 and 1 with an average switching rate jk- For 
a single fluctuator, the noise auto-correlation function is 



s(t)s(t') = cxp(-27|t - t'\). 
the power spectrum is given by 



(5) 



SrtnM =f- / s(0) S (t)e l ^dt 

^ J-oo 

1 4 7ff 2 
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As is well-known, an ensemble of fluctuators with I/7 
distribution of their switching rates gives rise to 1// noise 
power spectrum [Tij . 



With the criteria 



g k cos 6>^7 fc 



(7) 



we can put the fluctuators into two categories, the fast 
ones (<) and slow ones (>). K = M + N where M 
(N) is the number of slow (fast) fluctuators. The fast 
and slow fluctuators have qual itatively different effects 
on the qubit time evolution |20l — I2-4J ] . Fast is synonymous 
with weakly-coupled or Markovian, as can be seen from 
Eq. [3 The fast fluctuators can be treated with Rcdfield 
theory and they give rise to exponential decay of phase 
coherence. On the other hand, slow is synonymous to 
strongly-coupled or non-Markovian and Rcdfield theory 
cannot be applied. In general, for classical Markovian 
noise or Gaussian noise, the dephasing rates can be re- 
lated to the noise spectral density and filter functions 
[25l - l27j . A list of filter functions for common pulsing se- 
quences can be found in Table I of Ref. [23] ■ 

In this classical noise model, decoherence is a result 
of averaging the unitary time evolutions over all the 
possible noise sequencies Sk(t). The quasi- Hamiltonian 
method allows us to carry out this averaging analyti- 
cally and treat the fast and slow fluctuators on an equal 
footing [ITj, [2l|. The qubit dynamics is described by a 
transfer matrix acting on the qubit Bloch vector, i.e., 
n(t) = T(t)n(t = 0), while the transfer matrix is gener- 
ated by a non-Hermitian quasi-Hamiltonian. In the case 
of a single qubit interacting with a single fluctuator, the 



quasi-Hamiltonian has the form 



BnL, 



T3 



L 



where Tj are Pauli matrices associated with the fluctu- 
ator, and Li are the 5*0(3) generators associated with 
the qubit Bloch vector. Note the classical two-valued 
fluctuating field is mapped into a spin-1/2 particle in 
this formalism. The transfer matrix is given by T(t) 



(xf\exp(—iH q t)\if) where 



[l;l]A/2 cor- 



respond to unbiased fluctuator. Exact diagonalization 
of H q is possible only for 9 = while perturbation ex- 
pansion can be used in general to calculate the transfer 
matrix T(t). 

Signals from common experimental pulsing protocols, 
such as energy relaxation (ER), Hahn spin echo (SE) 
and free- induction (FID), can be calculated with the 
quasi-Hamiltonian method as well [2l[. For these puls- 
ing schemes, the qubit is initially in the ground state and 
the probability of the qubit being in the excited state is 
measured at time t. In the ER scheme, a single tt pulse 
is applied at the beginning of the measurement. In the 
FID scheme, two tt/2 pulses arc applied, one at the be- 
ginning and the other at the end. The SE scheme has 
the two 7r/2 pulses as in the FID scheme and another tt 
pulse in the middle of the time evolution, i.e., t/2. For 
our qubit-fluctuators model, the pulsed signals are given 
by 



n SE (t) ~ e -(r a +rs)t 



M 



1 + ^ £ 1>" sm (5m cos 9t) 
M 

n F ro(t) ^e" (r2+r3)t cos-B t JJ cos (g m cos 9t) 



(8) 
(9) 



A I 



1 + eim tan (g m cos 9t) 



m—l 



(10) 



where the relaxation and dephasing rates are given by 

27^ sin 2 9 



i 4 7 2 +B 2 ' 



r 2 =^ + r 

M 

F.3 = £ 7m 
m—l 

N 2 2 

n— 1 



(ii) 

(12) 
(13) 

(14) 



Here e lm = j m / (g m cos9) and e 2m = g m /B are the 
small parameters of the perturbation theory. 

The decoherence rates Tx = 1/Ti, T 2 = 1/T 2 and 
are caused by the fast fluctuators and the equations for 
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them are consistent with Rcdficld results [25[. In the 
case of a single fast fluctuator, the decoherence rates for 
the echo experiment can be directly connected to the 
noise power spectral density, i.e., 1/Ti = sin 2 #Srtn(-Bo)j 
1/T = cos 2 #S RT n(0) and 1/T 2 = 1/2T X + I/T4. T 3 is 
entirely due to the slow fluctuators. 



III. DETERMINATION OF MODEL 
PARAMETERS FROM EXPERIMENTAL DATA 

For purposes of data analysis, it is necessary for us 
to specify a not completely general but yet still flexible 
model for the noise. Let d (7) = J2k=i $ (l ~ 7 fc ) be the 
distribution of rates and take g to be independent of k, 
i.e., gk = g- If there is a range of couplings then g in 
the following formulas can be regarded as an appropriate 
average coupling. This equal-coupling-strength or single- 
coupling-strength assumption should not be a severe lim- 
itation of our model as long as the standard deviation in 
the distribution of g^'s is small relative to g itself, and to 
the width of the distribution of 7fc's. For the specific case 
of fluctuating magnetic moments producing flux noise, 
the model is appropriate if the moments are all on the 
surface of the superconducting loop. We will comment 
further on this below. 

We will assume a broad noise spectrum by taking 



d( 7 ) 



for 7 min < 7 < 7max 

otherwise. 



(15) 



Here 7 max and 7 m i n are the upper and lower cuts of the 
fluctuators' switching rates. For 1/f noise we must have 
7min > and 7 max < 00 in order that the energy den- 
sity of the noise be finite. The power-law assumption is 
often useful for analyzing experimental data, though the 
method used to solve the model itself is capable of treat- 
ing arbitrary distributions. Note s = gives 1/f noise. 
Under those assumptions, the pulsed signals are given 

by 



n ER {t) ~e- rit 
n SE (t) ^ e -0W3)* 



1 + — sin(7 c i) 

7c 



(16) 
(17) 



"FID 



(t) ~e-( r2+r ^cos.B icos M (7 c <) 



r 3 

1 H tan( 7c i) 

7c 



(18) 



where 7 C = g cos 9 is the critical coupling strength. The 
task of data analysis is then to determine the five intrin- 
sic parameters g, a, s, 7m i n , and 7 max from observations 
of the pulsed signals tier, "fid and ngE- The formulas 
show that the ER and SE signals at different working 
points alone are enough to fully determine all the five 
parameters, at least in principle. The FID data provide 
consistency checks and, crucially, to find the number of 
slow fluctuators at various working points. We have an- 



alyzed data from Ref. @, Q and all results are listed in 
Table H 

For the ease of analysis, it is convenient to define 



$SE (t,9) 



n SE (0) 



n SE {0 = tt/2) 

„-T 3 t 



1,^3. , 

1 H sin 7c r 

7c 



(19) 



where $se (t, 6 = 0) correspond to the ' pha se-memory 
functional' defined by other authors pil. |28|| . Note Fi 
has only weak dependence on the working point 8, thus 
r 2 (0) ~ T 1 (6 = tt/2)/2, and it drops out in Eq. [19] 
Similarly, we define $fid or the FID signal 



$fid {t, 1 



"fid (0) 



"fid (0 = tt/2) 
e -r 3 t 



1 H tan 7 c t 

7c 



cos M ( 7c t). (20) 



It is important to note that ^fid has explicit dependence 
on the number of slow RTN fluctuators M. 

We note the scaling parameter s has significant ef- 
fect on the working point dependence of the decoherence 
rates, especially for T^. In the case of 1/f noise, s = 
and we have 



ag 2 sin 2 9 _ x f 2j : 
1 1 ^ tan 



B 

r Tl 1 a 

r 2 ^ T + -7c 



Bn 



T 3 ~a(7 c - 7 m i„). 

Note r 3 has linear relationship to 7 C in this case. 
If s = 1, 



(21) 

(22) 
(23) 



log 



^+47c 2 

Tmax 



r 3 ~"( 7 2 - 7 2 . ). 

-j 2 ^ ' c ' mm/ 

In general, 

2asin 2 05 2 7 s+1 / s + 1 8 + 1 



(24) 

(25) 
(26) 



4T 2 

n 



1 



2(a 

(7c s+ 



Tj(7n 



7c 



s-11 



7 S+1 ). 

/ mm / 



(27) 
(28) 
(29) 



Here 2-F1 is the hypergeometric function. 

In flux qubit experiments, the working point is ex- 
perimentally tunable by varying the applied flux and 
7c = g cos 8 changes accordingly. Thus a plot of T 3 ver- 
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FIG. 1. (Color Online) Echo phase memory functional $sb 
data in Fig.4A of Ref.[3|. We fit the 36 data points (open 
circle) to Eq |19l (black solid line) and Gaussian model $f B = 
e~ F G t (blue dashed line) respectively. 



sus 7 C would unambiguously determine the distribution 
of the fluctuators d(-f). 



A. coupling constant g 

To extract the coupling constant g, we fit the phase 
memory functional to Eq. 1191 as seen in Fig. [TJ Thus for 
each working point 9, we can extract two numbers from 
the fitting, i.e. ^3(6) and 7c (0)- For example, the data 
in Fig. [JJwere taken at working point cos# = 0.18, and 
the corresponding T 3 = 0.99 MHz and 7 C = 2.1 MHz. 

In Ref. 0, H[ , the same data is fitted to a Gaussian 
noise model where the Gaussian flux fluctuation assumes 
a 1/f noise spectral density, i.e., S$(u>) = A/ui. In this 
case, the phase memory functional takes the Gaussian 
form $!? E = e~ r c* . As seen in Fig. [TJ both models 
fit the data well and it is is unclear which model is bet- 
ter. Non-Gaussian behavior manifests itself unambigu- 
ously with 'plateaus' in the phase memory functional $se 
[13, d^. The rise in the longer time in Fig. [JJcould be the 
onset of such 'plateaus'. A cleaner sample with fewer sur- 
face spins (smaller a) would help to make the 'plateaus' 
more visible, which would then distinguish the present 
model from the Gaussian model [2lj |. 

With data at different working point 6, we can plot -f c 
versus cos 9. The coupling constant g is the slope, as seen 
in Figj2ja) 0. Similar data analysis for Ref. has been 
carried out in Ref. 1211] - 



In the data analysis, we first extract Fq from Ref. 7, 8], then 

reproduce "^cfg from the Gaussian formula. Given fits the 

real experimental data well, we fit the reproduced to Eq. If 91 
to obtain and j c at different working points. 
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FIG. 2. Fitting of data from Ref.@]. Both 7c and T 3 are 
fitted from <&se at various working point, (a) Critical rate 
7c versus the working point cos 9. The slope is the coupling 
constant g. (b) Linear regression to Ts = q(7 c — 7min). 



B. noise intensity a, noise index s and lower cut 

7min 

The functional form of T3 allows us to determine a, s 
and 7 m in, as seen in Eq. [23 

Fitting the data from Ref. @ we get clean linear de- 
pendence of 7c, as seen in Fig. [UJb). A similar result has 
been obtained in Ref. 21 for the data from Ref. 0. 
This is a sign of 1/f noise in the environment (s = 0). 
7m in is the intercept of the linear fit. Unfortunately, its 



accuracy depends strongly on the quality of the data at 
low 7c, and the data are lacking in that region. Hence 
there is considerable uncertainty in the fitted value of 
7min- The noise intensity a can be retrieved from the 
slope, in the case of 1/f noise, as seen in Eq. [231 

If s is taken as a fitting parameter as well, we get s = 
-0.005 for Ref. @ and s = -0.07 for Ref. 0. Thus 
the data is consistent with 1 /f noise and we adopted the 
linear fit as in Eq. [231 When the noise power spectral 
density deviates farther from 1//, there are indications 
that as s decreases from 0, the pure dephasing time 
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B 


A/h (GHz) 


5.445 


3.9 


e/h (GHz) 


- 1 


- 1.1 


I P (nA) 


160 


370 


r (fjLTa) 


1 


2 


r[ ex) (MHz) 


0.65 


7.1 


g/h (MHz) 


9.5(0.2) 


68(2) 


7min (MHz) 


0.05(0.01) 


0.11(0.05) 


a 


0.77(0.006) 


0.754(0.003) 


rf h) (MHz) 


0.02(0.01) 


1.4(0.1) 


TO (/i B ) 


1.7(0.03) x 10 J 


5.2(0.2) x 10 a 


= lite) ($ /ite) 


3.5(0.2) x 10- 11 


3.3(0.2) x 10- 1 



TABLE I. Noise characteristics extracted from Ref.Q (col- 
umn A) and Ref.[^ (column B). The numbers in the top por- 
tion are experimental data while the ones in the lower portion 
are derived parameters from the model. The numbers in the 
parenthesis are standard deviations from linear regression. 



also decreases I2J 



C. upper cutoff 7 max 

For general s, T\ can be expressed in terms of hyperge- 
omctric function in s. In the case of 1/ / noise, s = and 
we have Eq. 1211 Since 7 ma x is the only unknown in the 
equation (ri is experimentally measurable and the other 
quantities can be derived from the experiment), 7 max can 
be extracted, at least in principle. But for the data in 
Ref. @, we are unable to back out 7, nax . One finds 
tan~ 1 (27 max /i?o) has to be greater than 7r/2 to validate 
the equation. 

The lack of experimental accuracy might not help for 
this self-inconsistency. One possible remedy is that there 
is some other source of high frequency noise, other than 
RTN, to cause relaxation. Thus the experimentally ob- 
served Ti, as denoted by T^ x ^ in Table U is actually 
greater than = irag 2 /2Bq which only includes the 
energy relaxation due to the fluctuators. Here is 



defined with j„ 



oo and 6 = ir/2 for Eq. |2T] 



D. number of fluctuators K 

Since most of the parameters in the model can be de- 
rived from the experimental data for tt-se and tier, we 
may use the data for npm to get a value for M(0), the 
number of slow fluctuators. As seen in Eq. the FID 
signal has explicit dependence on M. It is easiest to 
analyze this using the logarithm of the phase memory 
functional /C E (t). 

For the Echo signal, the logarithm of the phase memory 
functional JC E (t) can be expanded in terms of j c t <C 1 



and we have 

IC E (t) = 



log $SE 

r 3 7 c 2 i 3 /6, 

rjt 2 /2, 



ict > r 3 /7 c , 
j c t < r 3 /7 c . 



(30) 



Similarly, we define JC p for the envelope of FID signal. 
In the limit of j c t <C 1, 



IC F (t) 



log $FID, 



(31) 



Note at small times ( 7c i <C r3/ 7c ), both JC E (t) and 
lCp(t) are quadratic in time. If the waiting time for 
the first plateau is too long comparing to the damping 
time td — l/^s, both phase memory functionals will 
assume Gaussian shape. In the experiments 0, II], the 
authors fit to the quadratic terms in Eqs. 130] and |3"T1 
What they called T^ E and T^ F correspond to 



and V( A/ 7c 



r|)/2. 



In Ref. 0, ~ 8r^ £ . This linear dependence 

is expected as long as the qubit is not operated ex- 
tremely close to the optimal point 9 = tt/2, as can be 
seen from Eq. [23] Thus the number of slow fluctuators 
M ~ 64r|/7^ is of the order 10 for the working points 
in the experiment. To be more specific, 



M - 64 x ol 



38. 



(32) 



It should be evident that this is a rough estimate. How- 
ever, it is unlikely to be off by order of magnitude and 
we assert that M < 100. 



E. effective magnetic moment and power spectrum 

As a consistency check, we find the magnetic moment 
associated with the fluctuators and the total spectral den- 
sity. The change in flux due to spin in the SQUID loop 
is 



A$=^TO 



(33) 



where to is the effective magnetic moment of the spin, /xo 
is the magnetic constant and r is the radius of the loop. 

For flux qubit, we have 



21 J 



(34) 



where I p is persistent current along the qubit loop. Thus 

rg 



2/ioip 



(35) 
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The noise power spectrum density is 



S 1/f (uj) 



tan 



5 , RTN(^)rf(7)^7 

1 / 27mi 



tan 



27min 



' 2lu 



(36) 



With Eq. 1341 the noise power spectrum density in terms 
of flux is 

S*(w) = S 1/f (u)/4I$. (37) 
These derived results are listed in Table [H 

IV. CONCLUSION 



and find that the number of slow fluctuators is surpris- 
ingly small, less than 100. If we assume that the noise 
arises from magnetic clusters on the surface of the super- 
conducting loop, then the size of the magnetic moment of 
the clusters is quite large, of the order of 1000 to 5000/is. 

These results appear to be rather surprising. How- 
ever, a recent analy sis |30l of noise measurements on dc 
SQUID inductance [18[ suggests that the predominant 
noise sources are large magnetic clusters and that such 
clusters would give rise to 1/ f-like noise. This gives rise to 
a consistent picture of two quite different qubits analyzed 
in two quite different ways. 

The quasi-Hamiltonian method is applicable to more 
complex systems as well, such as interacting two-qubit 
systems [3l| . An extension of the current work would 
be to examine the more recent experiment where the de- 
phasing of two inductively coupled flux qubits are studied 



We have given a method for extracting the properties 
of the two-state fluctuators that cause decoherence in 
superconducting qubits. This method applies when the 
working point of the qubit can be varied, as is possible in 
flux qubit set-ups. The shape and strength of the noise 
spectrum can be determined from qubit measurements, 
and an estimate of the total number of active slow fluc- 
tuators can be obtained. We analyze two experiments 



ACKNOWLEDGMENTS 

We thank Robert McDcrmott, J. S. Tsai and K. 
Kechedzhi for helpful discussions and Zhigeng Gcng for 
teaching the authors to use R (the programming language 
and software environment for statistical computing and 
graphics). This work is supported by the DARPA/MTO 
QuEST program through a grant from AFOSR. 



[1] J. Clarke and F. Wilhelm, Nature 453, 1031 (2008) 

[2] Y. Makhlin, G. Schon, and A. Shnirman, 

|Rev. Mod. PnysT, 73, 357 (May 2001) [7] 

[3] Y. Nakamura, C. Chen, and J. Tsai, 

Physical Review Letters] 79, 2328 (Sep. 1997) V. Bouch- [8] 

iat, D. Vion, P. Joyez, D. Esteve, and M. H. De- 

voret, |Physica Scripta| T76, 165 (1998) Y. Nakamura, 

Y. Pashkin, and J. Tsai, Nature 398, 786 (1999) D. Vion, [9] 

a. Aassime, a. Cottet, P. Joyez, H. Pothier, C. Urbina, 

D. Esteve, and M. H. Devoret, [Science (New York, NY/)] [10] 

296, 886 (May 2002) 

[4] Y. Yu, S. Han, X. Ch u, S.-I. Chu, and Z. Wang, 

|Science (New York, NYQ] 296, 889 (May 2002) [11] 

J. Martinis, S. Nam, J. Aumentado, and C. Urbina, 

|Physical Review Letters] 89, 117901 (Aug. 2002) 

A. Berkley, H. Xu, R. Ramos, M. Gubrud, F. Strauch, 

P. Johnson, J. Anderson, A. Dragt, C. Lobb, and [12] 

F. Wellstood, ISciencel 300. 1548 (2003) 

[5] J. E. Mooij, T. P. Orland o, L. Lev itov, L. Tian, van der [13] 
Wal, C. H., and S. Lloyd. IScience 285, 1036 (Aug. 1999) 
C. H. van der Wal, A. C. J. ter Haar, F. K. Wil- [14] 
helm, R. N. Schouten, C. J. P. M. Harmans, T. P. Or- [15] 
lando, S. Lloyd, and J. E. Mooij, \ibid.\ 290, 773 (Oct. 
2000) J. Fried man, V . Patel, W. Chen, S. Tolpygo, 
and J. Lukens, [Nature] 406, 43 (Jul. 2000) I. Chiorescu, [16] 
Y. N akamura, C. J. P. M . Harmans, and J. E. Mooij, 
|Science (New York, N.Y)] 299, 1869 (Mar. 2003) [17] 

[6] H. Paik, D. Schuster, L. Bishop, G. Kirchmair, G. Cate- 
lani, a. Sears, B. Johnson, M. Reagor, L. Frunzio, 



L. Glazman, S. Girvin, M. Devoret, and R. Schoelkopf, 

Physical Review Letters] 107, 1 (Dec. 2011) 

F. Yoshihara, K. Harrabi, A. O. Niskanen, Y. Nakamura, 

and J. S. Tsai, |Phys. Rev. Lett.| 97, 167001 (Oct. 2006) 

K. Kakuyanagi, T. Meno, S. Saito, H. Nakano, 

K. Semba, H. Takayanagi, F. Deppe, and A. Shnirman, 

|Phys. Rev. LetE] 98, 47004 (Jan. 2007) 

F. C. Wellstood, C. Urbina, and J. Clarke, Appl. Phys. 

Lett. 50, 772 (1987) 

V. Foglietti, W. J. Gallagher, M. B. Ketchen, A. W. 
Kleinsasser, R. H. Koch, S. I. Raider, and R. L. Sand- 
strom, Appl. Phys. Lett. 49, 1393 (1986) 
R. C. Bialczak, R. McDermott, M. Ansmann, 
M. Hofheinz, N. Katz, E. Lucero, M. Neeley, A. D. 
O' Connell, H. Wang , A. N. Cleland, and J. M. Marti- 
nis, [PhysTRevTLett] 99, 187006 (Nov. 2007) 
R. H. Koch, D. P. DiVincenzo, and J. Clarke, 
|Phys. Rev. LetE] 98, 26700 3 (Jun. 2007) 
L. Faoro and L. B. Ioffe, |Phys. Rev. Lett.| 100, 227005 
(Jun. 2008) 

R. de Sousa, [Phys. Rev. B| 76, 245306 (Dec. 2007) 
S. Sendelbach, D. Hover, A. Kittel, M. Muck , J. M. Mar- 
tinis, and R. McDermott, |Phys. Rev. Lett. 1 100, 227006 
(Jun. 2008) 

F. Yoshihara, Y. Nakamura, and J. S. Tsai, 
Physical Review B| 81, 132502 (Apr. 2010) 
S. Gustavsson, J. Bylander, F. Yan, W. Oliver, F. Yoshi- 
hara, and Y. Nakamura, Physical Review B 84, 014525 
(Jul. 2011) 



7 



[18] S. Sendelbach, D . Hover, M. Muck, and R. McDermott, 

|Phys. Rev. LetE1 l03, 117001 (Sep. 2009) 
[19] S. Kogan, Electronic Noise and Fluctuations in Solids 

(Cambridge University Press, 2008) 
[20] R. Joynt, D. Zhou, and Q.-H. Wang, Int. J. Mod. B 25, 

2115 (2011) 

[21] D. Zhou and R. Joynt, |Phys. Rev.~A"l 81, 10103 (Jan. 

2010) 

[22] D. Zhou, A. Lang, and R. Joynt, [Gjuant. I nfo. Processing 

9, 727 (Mar. 2010) 
[23] E. Paladino, L. Faoro, G. Falci, and R. Fazio, 

|Phys. Rev. Lett .| 88, 228304 (May 2002) 
[24] Y. M. G alperin, B. L. Alt shuler, J. Bergli, and D. V. 

Shantsev, |Phys. Rev. LetE] 96, 97009 (Mar. 2006) 
[25] C. P. Slichter, Principles of Magnetic Resonance, third 

edit ed. (Springer, New York, 1996) 



[26] J. M. Martinis, S. Nam, J. Aumentado, K. M. Lang, and 
C. Urbina, |Phys. Rev. B| 67, 94510 (Mar. 2003) 

[27] L. Cy winski, R. M. Lutchyn, C. P. Nave, and S. Das 
Sarma, |Phys. Rev. B| 77, 174509 (May 2008) 

[28] Y. M. Galperin, B. L. Altshule r, J. Bergli, D. Shantsev, 
and V. Vinokur, [Phys. Rev. B| 76, 64531 (Aug. 2007) 

[29] S. Anton, C. Mueller, J. Birenbaum, S. O'Kelley, 
A. Fefferman, D. Golubev, G. Hilton, H. Cho, K. Ir- 
win, F. Wellstood, G. Schoen, A. Shnirman, and 
J. Clarke, Arxiv prepr int arXiv:1111.7272, 1(2011) 
larXiv:arXiv:1111.7272vTl 

[30] K. Kechedzhi, L. Faoro , and L. B. Ioffe(2011) 
|arXiv:1102.3445 [cond-mat]| 

[31] AT De, AT Lang, D7 Zhou, and R. Joynt, 
|Physical Review A| 83, 42331 (Apr. 2011) 



